Polaron-like effects in a one-dimensional optical lattice 
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We study a highly imbalanced Fermi gas in a one-dimensional optical lattice from the polaronic 
point of view. The time-evolving block decimationg algorithm is used to calculate the ground 
state and dynamics of the system. We find qualitatively similar polaronic behaviour as in the 
recent experiment by Schirotzek et al. [T] where radio-frequency spectroscopy was used to observe 
polarons in three-dimensional space. In the weakly interacting limit our exact results are in excellent 
agreement with a polaron ansatz, and in the strongly interacting limit the results match with an 
approximative solution of the Bethe ansatz, suggesting a crossover from a quasiparticle to a charge- 
density excitation regime. 

PACS numbers: 



I. INTRODUCTION 

Impurity problems are essential in determining low- 
temperature properties of condensed matter systems. A 
well-known example is an electron moving in a crystal lat- 
tice and interacting with its surrounding ions [SJ. This 
creates lattice polarization and deformation which is car- 
ried with the electron through the material. Interactions 
with the ions create an effective potential for the electron 
and try to slow it down, which can be modelled as an ef- 
fective mass for the electron. Another famous impurity 
problem is the Kondo effect [5] where the scatterings of 
conduction electrons with magnetic impurities give rise 
to electrical resistivity. Ultracoldic atom gases provide an 
excellent playground for the study of impurity problems 
due to the controllability of the system parameters. For 
instance, two-component gases can be realized by using 
two different hyperfme spin states of alkali atoms. The 
population imbalance can be controlled by transferring 
particles between the two hyperfme states with radio- 
frequency pulses. The inter-component interaction can 
be tuned by Feshbach resonances [3]. 

Recently, polarons were observed in an ultracold 
atomic Fermi gas [TJ |S]. In [JJ radio- frequency spec- 
troscopy was used to measure a sharp quasiparticle peak 
solely for the minority component in the highly imbal- 
anced gas of 6 Li atoms. Theoretically, various approaches 
such as Monte Carlo studies [5J [7] , T-matrix approaches 
[8HTU] and variational ansatze [5J [TTJ [T2J have been used 
to model the phenomenon. A variational ansatz by F. 
Chevy pj] has explained the qualitative features of the 
experiment. In |12j . variational ansatz describing also the 
BEC side of the Feshbach resonance, including molecular 
formation, was proposed, and the ground state properties 
matched well with the Monte Carlo studies [5J. 

In one-dimensional optical lattices, exact methods 
provide straightforward approaches for highly polarized 
gases [T3HTT] . In this article, we use an exact numeri- 
cal method to investigate highly imbalanced Fermi gases 
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from the polaronic point of view. We study both the 
ground state properties and the dynamics, and make a 
comparison to two different approximative solutions. In 
the weakly interacting limit the results are compared to 
the polaron ansatz of and in the strongly interact- 
ing limit the results are compared to an approximative 
solution of the Bethe ansatz. In both limits, an excel- 
lent agreement between the exact and approximative so- 
lutions is found. However, although we analyze numeri- 
cal results using a polaron ansatz, we do not necessarily 
claim the excistence of a polaron in a one-dimensional 
system: well-defined polaronic quasiparticles may not ex- 
ist due to the one-dimensional nature of the system. We 
discuss implications of our results to three-dimensional 
systems. 

In Section [IT] we consider ground state properties of 
highly imbalanced gases. In Section |Hl| we study radio- 
frequency response of the ground states. In Section [TV] 
limitations of our model as well as the connection to 
higher dimensional systems is discussed. In Section [V] 
we present the conclusions. 

II. EXACT SIMULATIONS 

For atoms in an optical lattice, the physics is well cap- 
tured by the Fcrmi-Hubbard model 

H = -J ^<y £ i° + U Y1 ™*T%> (!) 

where J is the hopping energy, U is the on-site interaction 
strength between different spin components a G {T,|}, 
Ci a (cJ CT ) annihilate (create) a fermion for a site i with 
spin a and hi a = c\ a ii . Here we consider only attrac- 
tive interactions U < 0. In recent experiments in [TJ [5] 
in highly polarized gases qualitative features have been 
explained with the variational ansatz 

I*> = ^o4J^>tI0H 

q<k] r ,k>k r F 
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FIG. 1: The polaron energy E p as a function of the in- 
teraction strength (—17). The weakly interacting regime is 
well described by the ansatz of Eq. |2| (red solid line) and 
the strongly interacting limit with the Bethe ansatz with 
—U/J —¥ oo (green dashed line). The number of sites is 
Nl = 40 and the atom numbers are Nf = 20, N± = 1. 



where \FS) refers to a filled Fermi sea, |0) to the vacuum 
state and k F is a Fermi momentum for spin f particles. 
In the first term the minority particle lies at the bot- 
tom of the band and majority atoms form the Fermi sea. 
In the second term the minority atom is scattered out 
from the lowest momentum state and the filled Fermi sea 
of majority particles is broken to particle-hole excita- 
tions. The variational coefficients ipo and ipkq are found 
by minimizing the energy. The energy difference between 
the interacting and non-interacting ground states is often 
named as the polaron energy E p = E g —E'g° n ~ mt because 
this energy difference corresponds to the energy needed 
to add a single impurity. Note that with our definition 
E p < for attractive interaction strengths. The quantity 
Z = \(fo\ 2 is a measure of the quasiparticle weight of the 
polaron. 

In order to see how well the ansatz describes the system 
we have calculated the exact ground state using the time- 
evolving block decimation (TEBD) algorithm 18J which 
allows us to determine several observables. The polaron 
energy is the expectation value E p — (g\H \g) — Eg° n ~ mt 
and the quasiparticle weight is the square of the inner- 
product Z = \{g\g non - int )\ 2 where \g) and \g n ° n ~ mt ) 
are ground states for the interaction strengths U < 
and U = 0, respectively, both with same particle num- 
bers. In Fig. [T] the polaron energy E p and in Fig. [2] 
the quasiparticle weight for the ansatz and exact numer- 
ics are shown. For small interaction strengths [/ < 5J 
the results are in good agreement but start to deviate for 
strong attractions. 

In three-dimensional free space the breakdown of the 
ansatz (pi) in the BEC limit has been shown experimen- 
tally [1] and described theoretically Q2]. In the limit 
l/(fcpa) <C 1 the emergence of molecular states is ex- 
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FIG. 2: The quasiparticle weight Z = \(po\ 2 as a function of 
the interaction strength — U. Parameters are as in Fig. [I] 



pected. However, in a one-dimensional optical lattice 
molecules are not formed even for strong, attractive in- 
teractions — U ^> 1. 

M. Punk et al. (H] provided an ansatz for the molec- 
ular regime. In their ansatz the dominant contribution 
comes from the term 



E M-A°to\ FS >^)u (3) 



where <pk are variational coefficients. The summation is 
now restricted to momenta above the Fermi momentum 
k^ F of majority particles. The momentum distributions 
for minority and majority components in our case are 
shown in Fig. [3] The momentum distribution for mi- 
nority component is extremely small above the Fermi 
momentum k F , and therefore the ansatz for molecular 
regime in highly polarized gases does not improve the 
ansatz ([2| in one-dimensional lattices. However we are 
able to explain our results in the strongly interacting 
limit quantitatively with the Bethe ansatz (BA). In the 
limit —U/J — > oo the Bethe ansatz solution can be ap- 
proximatively mapped to spinless fermions |19j . The en- 
ergy becomes (see Appendix) 



E=U-2J cos(fcj) 

(4) 



For AT-f = 20, ATj, = 1 and Nl = 40 the polaron energy 
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FIG. 3: Momentum distributions (fik, a ) = (c ktr Bka-) for the 
interaction strength U = — 10. The number of lattice sites 
is Nl = 40 and the atom numbers are Nf — 20, N± — 1. 
The pairing emerges between atoms close to the two Fermi 
surfaces, no pairing between k and — k above k\, is present. 

results from the calculation 

E p = E(U, N t ,NJ - E non - mt (U = 0,iV t ,iV+) 
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= U - 2J^cos(fcj) 

( \ 

- -2J^cos(fc 3 ) - 2Jcos(fci) 

= U + 2Jcos(fci) + 2 Jcos(fc 20 ) ~U + 2.07 J, 

where E(U, N t , ATj.) {E nm - int {U = 0, iVf.iVj.)) is the en- 
ergy for the system with interaction strength U (U — 0) 
and particle numbers N^-,N± (Nf,N±). The BA result is 
shown in Fig. [T] with the exact result. In the strongly 
interacting limit — U > 7 J the Bethe ansatz solution is 
in good agreement with the exact solution. 

III. RADIO-FREQUENCY SPECTROSCOPY 

In the experiment pQ, radio-frequency (rf-) spec- 
troscopy was performed to the minority and majority 
spin components. In rf-spectroscopy one of the spin com- 
ponents, j" or J,) is coupled to a third spin state, denoted 
by 3, which is not populated and sufficiently weakly inter- 
acting with the initial states. Theoretically, the rf-field 
is well described in the rotating wave approximation by 
the Hamiltonian 

H° r} {t) = <?» " V >^ + e l %c, CT ), (6) 
j 

where the sum is over the lattice sites, f2 is the coupling 
strength and 5 is the detuning of the rf-field from the 
(7 — 3 transfer frequency. 



FIG. 4: Minority and majority spectrum for U = — 10. The 
inset is a zoom of higher detunings where a quasiparticle peak 
emerges. The vertical line is the polaron energy E p calculated 
from the exact ground state. The number of spin components 
are N-f = 20, N± = 1 and the number of the lattice sites is 
Nl = 40. Spectra for other interactions are qualitatively 
similar. 

From the variational ansatz ^ we have straight access 
to the rf-spectra through the Fermi golden rule 

I^J2\(f\^f\9)\ 2 S {ld HS-E f + E g ) 

f (7) 

where the summation is over all states with energies 
Ef, \g) is the ground state and E g its energy jTl. The 
first term gives rise to a narrow peak at the polaron en- 
ergy and behind that is a broad tail resulting from the 
term T mc (ipkq, 5) describing the incoherent part. Quali- 
tatively the measured spectra [1] matched with the varia- 
tional ansatz for the attractive interaction strenghts and 
close to the unitary limit. However, in the BEC side 
I/(fci?a) <C 1 the minority and the majority spectra over- 
lapped completely which signaled molecular pairing. The 
TEBD-algorithm allows us to calculate the full time- 
evolution of the system, and therefore we can evaluate 
the relation between the rf-spectra and the polaron-like 
state. To obtain the spectra, we first calculate the ground 
state of the Fermi-Hubbard Hamiltonian ([T]) and then the 
dynamics by operating on the ground state with the time- 
evolution operator exp(-iHt) consisting of the Hubbard 
Hamiltonian and the rf-field: 

H = H + H? f . (8) 

The spectra for the minority and majority components 
calculated using exact numerics are shown in Fig. [4] In 
the majority spectrum the main contribution is at zero 
detuning because of the unpaired atoms, but some atoms 
are transferred also when 5 ~ 8— 14 J. The minority spec- 
tra is highly peaked at S = 8.4 J, and the peak position is 
exactly at the polaron energy E p which is calculated from 
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FIG. 5: Minority component spectra for various imbalances. 
The interaction strength is U = —10 J. 

the exact ground state, and is very well approximated by 
Eq. ([5]). At larger detunings a long tail emerges but is 
cut due to the restriction of momenta in an optical lattice 

tO k G [— 7T, 7r]. 

In Eq. Q the particle-hole excitations give rise to the 
incoherent part of the spectra 

T mc ( Vkqi 5) = \^ kq \ 2 6M(6 - e q - k -e k + e q + e - \E p \), 

(9) 

where e k — — 2Jcos(/c) is the dispersion relation for non- 
interacting particles in a lattice. The minimum detuning 
which contributes to the minority spectrum is for q = k 
and gives 8 m i n — E p . The maximum detuning arises 
when the system has the hole in the bottom of the Fermi 
sea q = and the excitation lies at the van Hove singular- 
ity p — 7r i.e. S max = e-jT + - e + E p - e = E p + 8J. 
Now, the width of the spectrum is S max — S m i n = 8J. 
This is in good agreement with the exact spectra. 

In Fig. [5] we have varied the polarization. The spec- 
tra look similar for large spin imbalance N^/N^ < 0.2 
but for small spin imbalance the spectra get broadened 
and are shifted to lower detunings. Let us analyze the 
peak positions with the Bethe ansatz. The Fermi golden 
rule states that the peak position comes from the energy 
difference between the ground state and final state. For 
ground state with = 20 and N± = 2 the final state for 
the Fermi golden rule analysis has the particle numbers 
iV-f = 20, iVj. = 1 and N3 = 1. Therefore, the energy 
difference becomes 

AE=E G - E F 
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= 2U- 2J^cos(fc l )- 

( " \ (10) 

U - 2J^ COs(fci) - 2, J COs(k f ina i) 
= U + 2JC0S(/Cig) + 2JcOs(kf ina i). 

The lowest contribution to the spectrum arises when 



k final = ki which implies 

\E P \ - \AE\ = 2Jcos(fci 9 ) - 2Jcos(fc 20 ) « 0.15. (II) 

The energy difference \AE\ is smaller than the polaron 
energy \E p \ (energy difference in the case of a single mi- 
nority component), and therefore the spectrum shifts to 
lower detunings when polarization is increased. Further- 
more, we can make quantitative comparison of this Bethc 
ansatz result to the exact numerics. The distance be- 
tween the two peaks in the rf-spectrum for the cases 
= 1, 2 is around 0.1, which is in good agreement with 
the Bethe ansatz result \E P \ — |Ai?| = 0.15. 

IV. DISCUSSION 

We have shown that in a highly imbalanced Fermi gas, 
the ground state exhibits polaronic-like behaviour. The 
quasiparticle nature of this excitation, which can be de- 
duced from the non-zero quasiparticle weight of polaron, 
is confirmed by the rf-spectroscopy analysis. Two con- 
siderations are, however, in order. The first concerns 
the dynamical properties of the quasiparticle. From our 
analysis it is not possible to prove the stability of the 
polaronic-like particle propagation through the cloud of 
majority atoms. The polaron weight might split up into 
particle-hole excitations. The rf-spectroscopy, measuring 
the single particle spectral function, does not describe 
collective properties. A nonzero value for quasiparticle 
residue would correspond to dynamical stability of the 
quasiparticle: 

Urn |G+(k = 0,f)|^0. (12) 

The second consideration relates to the quasiparticle de- 
scription of the polaron. On general grounds it is well 
established that one-dimensional systems exhibit a "col- 
lectivization" of the excitation, the typical example be- 
ing the spin-charge separation in the balanced Fermi- 
Hubbard Hamiltonian. However it has been shown that 
in the case of an imbalanced gas the spin-charge sepa- 
ration is violated O Q31 OH HI]- We expect that the 
validity of the ansatz given in Eq. ([2|, and hence the 
validity of the polaron quasiparticle description even in 
ID case, in the limit of weak interaction, can be pictured 
as the extreme limit of the violation of the spin-charge 
separation. 

V. CONCLUSIONS 

We have investigated the ground state properties and 
the rf-spectrum of a highly imbalanced Fermi gas in a 
onc-dimcnsional optical lattice. The exact numerical re- 
sults can be explained in terms of the variational ansatz 
given in Eq. ^ and of Bethe ansatz equation (in the 
limit of — U I J — > 00) for —U < 5 J and — U > 7 J, respec- 
tively. Our results suggest the possibility of the existence 
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of a polaronic quasiparticle, further analysis is however 
required to investigate its dynamical properties. More- 
over, we would like to point out that the setup proposed 
here is well within reach of the current experimental ca- 
pabilities. Our analysis provides exact numerical results 
for comparison with future experiments, as well as effec- 
tive physical interpretation for both the weak and strong 
interaction limits. 
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Appendix 

Eq. Q can be obtained considering the BA solution 
for the open-boundary conditions (OBC) Fermi-Hubbard 
model in the limit U/J — > oo. Analogously to the calcu- 
lations performed in [115], it is possible to prove that the 
excitations of the system, for a repulsive interaction, can 
be described in terms of N = N-f + N± spinlcss fermions 
with energy and momenta given respectively by 

N 



E = -2jV]cosfcj 



— Ij IjeN,j=[l...N]. 



(13) 



The expression for kj in Eq. ( 13 ) is derived from the BA 
equations for [5^. Nf up, N± down electrons on Nl sites 
BA equations can be written as 



2Lkj =2wlj 



2k 



E 

/3=1 



$ 2 



,sin(fcj) - \ p 



$ 2 



sin(fcj) + A/3 



(14) 



and 



$ 2 



E 

=2tt J a + 



t X a - sin(fcj) 



$ 2 



t \ a + sin(fcj) 



u 

A« + A/j 
u 



(15) 



where j = 1, . . . , N t + N, , a = 1, . . . , N h Ij , J a € N, 
U/J = u, $(.t) = 2 tan -1 (2 a;) and X a are the spin ve- 
locities. 

The distribution of Ij should correspond to a condition 
where the energy is minimized. For a balanced gas at 



half filling, the energy minimization condition is given 
by Ij = [1...L], leading to E = — 2 J J2j=i cos kj = 0) 

If we take into account the mapping U — > — U and 
that, in the limit U / J — > oo, the total number of pairs is 
equal to N±, the single-site basis states can be mapped 
according to the following scheme 



m>+*u>, I0)^it), 

leading to N = N L - (N± - iV t ), and hence 

Nt,-{Ni,-Nf) 



E= -23 



E 



COS kj 



(16) 



(17) 



which, taking into account the fact that J3j=i cos = 0> 
can be written as 



E = — 2 J cos kj 

3=1 



(18) 
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